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ABSTRACT 

A deep survey of the European Large Area ISO Survey (ELAIS) field in the southern 
celestial hemisphere (hereinafter SI) has been carried out with the Australia Telescope 
Compact Array (ATCA) at 1.4 GHz. The SI region, covering about 4 square degrees, 
has been also surveyed in the mid- and far-infrared (5-200 /im) with the Infrared Space 
Observatory (ISO). The radio survey provides uniform coverage of the entire SI region, 
with a sensitivity (5a) of 0.4 mJy over the whole area and 0.2 mJy in the center. To 
this sensitivity virtually all the radio counterparts of the far-infrared extragalactic 
ISO sources should be detected. This paper presents a radio sample-complete at the 
5a level-consisting of 581 sources detected at 1.4 GHz. Of these radio sources, 349 
have peak flux density in the range 0.2-1 mJy, forming a new homogeneous sample 
of sub-mJy radio sources. Due to its size, depth and multi-waveband coverage, the 
sample will allow us to study in greater detai l the sub-mJy radio sour ce population. 



The full catalogue will be available from bttp:/ / athena.ph.ac.uk/ 
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1 INTRODUCTION 

The European Large Area ISO Survey (ELAIS, Oliver et al. 
1997; Oliver et al. 1998 in preparation) is a collaboration 
between 20 European institutes. It involves a deep, wide- 
angle survey at high galactic latitudes, at wavelengths of 
6.7 /im, 15 /im, 90 |im and 175 /im with the Infrared Space 
Observatory (ISO). The 6.7 /jm and 15 /im surveys were 
carried out with the ISO-CAM camera (Cesarsky et al. 1997) 
to 5<t sensitivities of ~0.6 and 2 mJy, respectively. The 90 
fjxn and 175 /im surveys were carried out with the ISO- 
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PHOT photometer (Lemke et al. 1994), to a 5a sensitivity 
of ~60-80 mJy at 90 /im. With a sensitivity of ~60 mJy 
at 90 /im, ELAIS will be the deepest far-infrared survey 
performed with ISO. The survey has detected objects 5-10 
times fainter than IRAS in the 50-100 /im range, and 20-50 
times fainter than IRAS in the 10-20 /im range. The survey 
is divided into 4 fields (one of which, SI, is in the southern 
hemisphere), and covers a total area of ~ 13 deg 2 at 15 /im 
and 90 /im, ~7 deg 2 at 6.7 /im and ~3 deg 2 at 175 /im. 

ELAIS will allow us to study dust emission in normal 
galaxies and the evolution of star formation to high red- 
shifts, testing competing scenarios for the formation of el- 
liptical galaxies, exploring the IRAS galaxy population to 
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higher z, and possibly unveiling new classes of objects or 
unexpected phenomena. Since dust plays an important role 
in most of these goals, and since many of the thousands of 
galaxies detected in the ELAIS survey will be at high red- 
shift and probably obscured in the optical bands, radio ob- 
servations will play a crucial role in assessing the reliability 
of the ELAIS source list and in facilitating source identifica- 
tion in faint or empty optical fields. Additionally, radio data 
will be important in the optical identification phase because 
the spatial resolution of ISO is insufficient to identify un- 
ambiguously many of the faintest ISO sources. Even at 15 

the survey resolution is ~10 arcsec, while at 90 fira it 
is almost 1 arcmin, so there are multiple optical candidates 
within each error box. 

The three northern ELAIS fields (Nl, N2 and N3) have 
been surveyed in the radio band (at 1.4 GHz) with the Very 
Large Array (VLA) down to an average flux density limit 
of 0.25 mjy (5<r) over a total area of 4.22 deg 2 . Details of 
the VLA observations and a description of the catalogue 
containing 867 radio sources can be found in Ciliegi et al. 
(1998). 

In this paper we describe radio observations of the 
ELAIS field located in the southern celestial hemisphere, 
SI. This field, centered at a(2000) = 00 h 34 m 44.4 s , 5(2000) 
= -43° 28' 12", covers an area of the sky of about 2° x 2°. 
We observed the whole area at 1.4 GHz with the Australia 
Telescope Compact Array (ATCA) in the 6-km configura- 
tion (maximum baseline length), with a resolution of 8 x 15 
arcsec. We obtained uniform radio coverage over the whole 
SI region, with sufficient sensitivity (la ~ 80 /^Jy) to detect 
virtually all radio counterparts of the ISO galaxies. 

By combining these radio observations with the avail- 
able ISO data, we will investigate the radio/far- infrared cor- 
relation (Helou, Soifer & Rowan-Robinson 1985) in star- 
forming galaxies to a flux density significantly deeper than 
that reached by IRAS. When spectral informations become 
available for the optical counterparts of the ISO and ATCA 
sources in the field, we will investigate the trivariate IR- 
radio-optical luminosity function and its evolution. This 
study will be of fundamental importance for interpreting 
the source counts at different wave-bands and for elucidat- 
ing contributions from different classes of objects. Moreover, 
we will be able to determine the influence of obscuration on 
the inference of the star-formation rate, and to constrain the 
star formation history of the Universe to z ~ 1.0-1.5 with 
two complementary samples, selected in the infrared and ra- 
dio, within the same volume of the Universe (see Cram et 
al. 1998; Oliver, Gruppioni & Serjeant 1998). 

The radio sample is also important in its own right, 
since it constitutes one of the largest homogeneous sam- 
ples of sub-mjy radio sources, whose nature is still uncer- 
tain. Previous work has suggested that the sub-mjy radio 
source population is composed largely of star-burst galax- 
ies at moderate redshifts (i.e. Windhorst et al. 1985; Benn 
et al. 1993). However, there is little information about the 
nature and the true redshift distribution of the sub-mjy 
population, because only a small fraction of them have been 
optically identified and have measured redshifts (Benn et al. 
1993; Rowan- Robinson et al. 1993). Although it is commonly 
believed that most of the sub-mjy radio sources are associ- 
ated with star-forming galaxies, recent extensions to fainter 
optical magnitudes (and higher identification percentages) 



hint strongly at an increasing fraction of early-type galaxies 
(Gruppioni, Mignoli & Zamorani 1998). This suggestion, to- 
gether with apparently discordant results from optical iden- 
tifications of even fainter pjy radio samples ( Windhorst et 
al. 1995; Hammer et al. 1995; Richards et al. 1998) com- 
plicate the picture. The ATCA survey in the ELAIS field 
SI, together with the VLA surveys of the northern fields 
Nl, N2 and N3 (Ciliegi et al. 1998) will shed new light on 
the sub-mjy population and its connection to the faint blue 
and infrared star-burst galaxies. The sensitivity, areal cov- 
erage and multi-wavelength character (optical, near-, mid- 
and far-infrared and radio) of the ELAIS project will en- 
sure that it makes an important contribution to all of these 
problems. 

In section 2 of the paper we describe the radio obser- 
vations. Section 3 discusses the data reduction strategy and 
section 4 presents the source catalogue. In section 5 we dis- 
cuss the radio source counts and present our conclusions in 
section 6. 



2 RADIO OBSERVATIONS 

2.1 Observing Strategy 

By choosing an observing frequency of 1.4 GHz, we ensure a 
large instantaneous field of view (33 arcmin primary beam) 
and hence shorten the observing time necessary to cover the 
entire SI area. The sensitivity is also higher for the same 
integration time with respect to higher available frequen- 
cies (2.4, 4.86 and 8.44 GHz). Moreover, a frequency of 1.4 
GHz reveals well any recent star-formation activity, seen at 
this frequency as synchrotron radiation excited by supernova 
remnants (Condon 1992). 

We observed in continuum mode with a bandwidth B 
of 128-MHz. The sensitivity AS (mjy) is 

AS = 0.07555^ (NtB)~ 1/2 (1) 

where S 3ys (in Jy) is the system sensibility (350 Jy at 1.4 
GHz), t is the integration time (in minutes), B is the band- 
width (in MHz) and N is the number of baselines. 

2.2 The Mosaic Technique 

The ATCA is an east-west array which uses the earth's ro- 
tation to sample visibilities over elliptical loci in the spa- 
tial frequency (u, v) plane. Its control system allows imag- 
ing over a field or a source whose extent is larger than the 
primary beam by cycling through a grid of pointings on 
the sky, recording the visibilities for each pointing intermit- 
tently. This mosaic observing mode is an efficient way to 
obtain uniform quality imaging over a large field, and has 
been used to construct our survey. 

The grid of pointings in the mosaic was designed to yield 
a homogeneous radio source sample. We aim to have a de- 
tection threshold independent of source position, and there- 
fore require uniform noise over the area of interest. With 
the mosaic technique, images obtained with single pointings 
are combined together into a large image (mosaic) of the en- 
tire observed region. We used 'linear combination' mosaicing 
rather than the 'joint deconvolution' approach (described by 
Sault, Staveley-Smith and Brouw 1996) because the latter 
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Figure 1. The sky position and orientation of the ISO SI survey 
region. The outer rectangle is 2° X 2°. The filled triangles show 
the ISO-CAM 15 pointings, while the filled squares show the 
ISO-PHOT 90 fim pointings. The circles (drawn with a radius of 
20 arcmin) show the ATCA mosaic pointings. 

requires greater computing resources and is less suited to 
images dominated by compact sources. Linear combination 
mosaicing takes a pixel-by-pixel weighted mean of the sin- 
gle pointings, where the weights are determined from the 
primary beam response and the noise level in that pixel 
(Sault & Killeen 1995). The value of each pixel in the fi- 
nal mosaic depends on the pointing configuration and, in 
particular, on the grid spacing. For our survey, an optimal 
trade-off between uniformity of sensitivity and efficiency of 
telescope usage was obtained by adopting a rectangular grid 
of 7 x 7 fields, with a regular spacing between fields equal 
to 9fwhp arcmin (where Qfwhp is the full width 

at half power of the primary beam, 33 arcmin at 1.4 GHz). 
Figure 1 shows the sky position and orientation of the SI 
region. Circles drawn with a radius of 20 arcmin show the 
49 ATCA pointings. 

2.3 Observations 

The ATCA observations of the ELAIS field SI were per- 
formed on the nights of 1997 June 23, 24, 26, 27 and 28, 
in runs of 12 hours per night. The theoretical sensitivity is 
80 /iJy (la) over the whole area. Since the far- infrared and 
radio luminosities of normal galaxies are tightly correlated 
(Helou, Soifer & Rowan-Robinson 1985; Condon 1992), flux- 
limited samples of normal galaxies selected at far-infrared 
and radio wavelengths at the corresponding sensitivity are 
nearly identical. Given a 5a sensitivity of 60 mjy for ISO ob- 
servations at 90 /im, the far-infrared/radio correlation pre- 
dicts that 1.4 GHz radio counterparts of the faintest ELAIS 
sources should be detected at a flux density a little below 
0.5 mjy. Thus, the flux density limit reached by our obser- 
vations, corresponding to 5a ~ 0.4 mjy, should allow us to 



detect essentially all the ELAIS ISO-PHOT and most of the 
ISO-CAM galaxies. 

To optimize (u, v) coverage we organised the observa- 
tions by dividing the 49 fields in four blocks of 13 fields each 
and observing each block for 12h + 3h (one night plus one 
quarter of the last night). During this time the 13 fields were 
observed with dwell times of 20 seconds each, adding 3 min- 
utes for observing the secondary calibrator every hour (i.e. 
every 15 cycles of 13 fields). The field corresponding to the 
center of the SI area was visited in each block, to obtain 
deeper sensitivity in the central ~ 20 arcmin of the ELAIS 
field. 

The primary flux density calibrator was PKS B1934- 
638, whose flux densities at different frequencies are incor- 
porated directly in the calibration software. The primary 
calibrator was observed for 10 min at the beginning of each 
12h observing run. The source PKS B0022-423 was used as 
a phase and secondary amplitude calibrator. 

Before observing we checked for the presence of strong 
radio sources in the field, which could eventually compro- 
mise the sensitivity of the survey. Although the ELAIS SI 
field was not selected with radio observations in mind, for- 
tunately there are only a few bright radio sources within 
it, so that the target noise level was obtained without seri- 
ous problems of confusion. The Parkes-MIT-NRAO (PMN) 
Southern Survey (Griffith & Wright 1993) reveals one bright 
radio source (PMN J0042-4413) just outside the field, with 
a flux density of 1.2 Jy at 5 GHz. There are also one 0.4 Jy 
and one 0.2 Jy PMN source in the field, and several with 
lower flux densities. We have observed the whole SI area, 
expecting that not more than a small part of the mosaic 
might be degraded by the sidelobes of strong sources in or 
just outside the field. 

We used both receivers at 1.4 GHz, each with a band- 
pass of 128 MHz. Using the ATCA in this mode makes it 
hard to avoid strong interference in some channels for at 
least part of the observation, but the doubled bandwidth 
improves the sensitivity by a factor of v2- The optimal po- 
sitioning of the bands was obtained by centering one fre- 
quency at 1.344 GHz (1.280-1.408 GHz) and the other at 
1.472 GHz (1.408-1.536 GHz). 



3 DATA REDUCTION 

The data were analyzed with the software package 
miriad (Multi-channel Image Reconstruction, Image Anal- 
ysis and Display), which is standard software for radio- 
interferometric data reduction specifically adapted to the 
ATCA. 

In particular, miriad provides multi-frequency synthe- 
sis (mfs) algorithms (Sault & Wieringa 1994), which give 
the opportunity of producing images with improved (u, v) 
coverage by combining accurately the visibilities of individ- 
ual channels. Multi-frequency synthesis may be used with 
the ATCA because, in continuum mode, the correlator pro- 
vides a bandwidth of 128 MHz subdivided into 32 frequency 
channels (of 4 MHz each). The division of the wide pass- 
band into sub-channels reduces the effects of bandwidth 
smearing. Moreover, mfs fills in the (w, v) plane despite the 
modest number of ATCA antennas, since different observing 
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Figure 2. Grey scale image of the total 1.4 GHz mosaic in the ELAIS SI field. The square region corresponds to the area observed by 
ISO. 



frequencies produce different spatial frequency intervals Am 
and Av. 

Before calibrating the data, we flagged all baselines and 
fringe visibilities affected by any problems encountered dur- 
ing the observation (such as antenna slewing intervals, inter- 
ference, etc) using the tasks blflag, tvflag and uvflag. 
Primary calibrator data were flagged before calibrating, and 
the secondary calibrator after application of the band-pass 
and instrumental polarization calibrations obtained from the 
primary calibrator. Finally, the fully calibrated single field 
data were flagged. Each bandpass (1.344 and 1.472 GHz) was 
calibrated separately, according to the standard calibration 
technique (Sault and Killeen 1995). 

The complete imaging procedure (described in detail 
by Prandoni 1997) comprises imaging itself (performed with 
the miriad task invert), and several cycles of cleaning (task 
clean) and self-calibration (task selfcal) . The whole pro- 
cedure was applied separately to each frequency to account 
for spectral variations of each source, since flux density dif- 



ferences due to the intrinsic spectral slope could be signif- 
icant over the frequency interval of 2 x 128 MHz. For each 
field we produced individual images with the same geome- 
try and the same celestial coordinate as reference point, to 
be combined together into a single mosaic at the end of the 
reduction phase. 

The images are 1200 x 1200 pixels square, and for each 
field a 2400 x 2400 pixels square beam was produced. An 
average synthesized beam of 15 x 8 arcsec 2 was created for 
all the images. Cell sizes of 2.5 arcsec per pixel were used, to 
obtain good re-sampling of the (it, v) plane (at least 3 pix- 
els within each dimension of the synthesized beam). Natural 
weighting was used, to increase the sensitivity at the price of 
a decrease in spatial resolution. The images for each point- 
ing center were then CLEANed. A few preliminary deconvo- 
lution iterations (performed to find the model components 
for the self-calibration) were followed by a self-calibration 
phase and then a deeper cleaning phase. This cycle was re- 
peated twice. After subtracting from the visibility file the 
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Figure 3. The rms noise as a function of distance from the field 
center, measured along the two diagonals of the square defined 
by the SI field. Top: rms noise along the diagonal degraded by 
the sidelobes of PMNJ0042-4413; Bottom: rms noise along the 
diagonal not affected by the presence of the strong radio source. 
The solid lines (corresponding to lcr rms noise values of 80, 60 
and 40 fijj for r > 20 arcmin, 10 < r < 20 arcmin, and r < 10 
arcmin, respectively) represent the limiting fluxes assumed for the 
source extraction. 



components found by the clean algorithm, we flagged any 
residual visibilities still affected by interference. Then we 
performed the deconvolution on the new visibility data-set 
(after having recombined the previously subtracted compo- 
nents). Self-calibration was used to make additional correc- 
tions to the antenna gains as a function of time, to improve 
the image quality and to account for the fact that interpo- 
lated calibrator gains do not determine the antenna gains 
perfectly at each time step. 

To minimise the side-lobes of bright off- field sources, 
we performed a few iterations of the cleaning algorithm on 
a double-sized (four times in area) low resolution image of 
the field, subtracting the components found in the external 
parts of the image from the visibility file. The final images 
were obtained by adding the components found by the de- 
convolution algorithm to the residual image. This was done 
with the task RESTOR. Once good quality images of each 
field were obtained, we constructed the final mosaics using 
the linear combination task linmos (which also corrects for 
primary beam attenuation). Two mosaics were first created 
for each of the two frequency bands and then added, result- 
ing in a square 2° x 2° mosaic with almost uniform sensitivity 
over the full region. Figure ^| exhibits a grey scale image of 
the total mosaic of the entire SI area. 



Figure 4. Regions of different average rms noise defined for 
source extraction. The radii of the circles are 10, 20 and 50 ar- 
cmin for the central inner, the central outer and corner circles 
respectively, and the corresponding rms noise is 40 fijj, 60 /ijy 
and 100 fijy. The rest of the field has a constant rms of 80 fijy. 
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Figure 5. Areal coverage of SI field represented by the solid angle 
over which a source with peak flux S can be detected. 



3.1 Noise distribution in the mosaic 

The amplitude distribution of pixel values in the final mosaic 
corresponds to an almost Gaussian noise core plus a posi- 
tive tail due to actual sources. The standard deviation of 
the Gaussian which fits the distribution is ~80 /iJy, almost 
coincident with the theoretical rms value expected in the 
image. However, as is clearly visible in figure 0, the bright, 
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Figure 6. Distribution of peak flux densities for radio sources of 
the SI 1.4 GHz complete sample. 

extended source PMNJ0042-4413 located just outside the 
field of interest was strong enough (above 1 Jy) to produce 
sidelobes and hence limit the dynamic range over a part of 
the mosaic. 

To deal with this problem, we analyzed the noise prop- 
erties in the image as a function of radial distance alongtwo 
diagonals of the square defined by the SI field. In figure H we 
plot the variations of the rms noise with distance from the 
field center along the diagonal not affected by the sidelobes 
of PMNJ0042-4413 (top panel) and along the one degraded 
by them (bottom panel). Except in the corner disturbed by 
PMNJ0042-4413, where the rms noise is obviously higher 
than expected, the rms noise in the mosaic is uniform over 
most of the image. It is smaller in the center where a longer 
integration has been obtained, and almost flat over most 
of the survey area. Using figure ^| we defined regions with 
different but uniform levels of rms noise to characterise the 
sensitivity for source extraction. We adopted a la sensitivity 
of 40 fiiy in a circular area of 10 arcmin radius in the field 
center; 60 fiJy in a 10 < r < 20 arcmin annulus; 100 fiJy in 
the area delimited by a 50 arcmin radius circle centered on 
the noisy corner and two sides of SI; and 80 piJy over the 
rest of the field. Figure [l| illustrates the location and extent 
of these regions. 

Using these sensitivities we obtained the integral distri- 
bution of the rms noise in the image and the detectability 
area as a function of flux density. In figure ^ the solid angle 
Q(S) over which a source with a peak flux density S can be 
detected is plotted as a function of flux density. The 'step' 
structure of the areal coverage is due to the fact that the 
field is divided in separate regions of constant rms noise. 



4 THE SOURCE CATALOGUE 
4.1 Source detection 

The criterion for including a source in the catalogue is that 
its peak flux density is > 5 times the average rms value 



in its region of the image. Source extraction was performed 
using the MIRIAD task imsad (Image Search and Destroy), 
which searches for islands of pixels above a given cutoff and 
attempts to fit Gaussian components to the islands. Source 
parameters, derived by least-squares bi-dimensional Gaus- 
sian fitting, are the right ascension and declination of the 
island centroid, peak flux density, integrated flux density, 
deconvolved (from the beam) major axis full width at half 
maximum FWHM (arcsec) , deconvolved minor axis FWHM 
(arcsec) and deconvolved position angle (degrees). 

As discussed by Condon (1997), the results of Gaus- 
sian fitting can be unreliable for sources with low signal-to- 
noise ratios. Thus, we used imsad to extract all the sources 
whose peak flux, S p , was greater than 4 times the local rms 
value and then, for sources with 4a < S p < 7a, we derived 
the peak flux density by second degree interpolation (task 
maxfit) and the total flux density by integrating the image 
value in a rectangle around the source. Only sources with 
a maxfit peak flux density > 5a were included in the final 
sample. For the other parameters (major axis, minor axis 
and position angle) we retained the values given by IMSAD. 

4.2 The catalogue 

Our method selected 581 sources with peak flux density > 
5a, over an area of 4 deg 2 . Of these, 31 appear to have two 
or more components. Figure ^| illustrates the distribution of 
peak flux density for the 581 sources in the catalogue. This 
distribution shows that the majority of the radio sources 
in our sample (~ 60%, 349/581) have peak flux density in 
the sub-mjy region. Thus, our data-set provides a large and 
statistically significant sub-sample of radio sources fainter 
than 1 mJy. 

We have considered as potential doubles all the sources 
separated by less than about twice the FWHM of the synthe- 
sized beam (~ 30 arcsec) and having approximately equal 
flux densities in the two components (ratio < 2.5). For triple 
or multiple sources the distance between each component 
and the probable nucleus has been considered, while the 
flux density ratios have been computed between the various 
components excluding the nucleus. All the components sat- 
isfying the adopted criteria have been assumed to form a 
unique, multiple source. Figure Q shows contour images of 
the sources classified as double or multiple in the catalogue. 

The catalogue contains a total of 621 components, 
and reports the source name, the peak flux density S p 
(in mJy), the integrated flux density St (in mJy), the 
source position (right ascension and declination at equinox 
J2000), the FWHM of the major and minor axes, the po- 
sitional angle of the major axis, the signal-to-noise ra- 
tio and a character as a comment about the deconvolu- 
tion outcome (D=deconvolution of source size from beam 
OK, P=deconvolution gave result close to point source, 
F=deconvolution failed, E=multiple source). Whenever the 
deconvolution of the source size from the beam failed (F), 
the fitted parameters (instead of the deconvolved ones) are 
reported in the catalogue. For double or multiple sources 
the components are labeled 'A', 'B', etc., followed by a 
line labeled 'T' in which parameters for the total source 
are given. The position of the total source has been com- 
puted as the flux-weighted average position for all the com- 
ponents. In Table 1 the first page of the catalogue is shown 
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Table 1: Tb 



Name S» 
(mjy) 

ELAISR20J002830-424456A 7.182 

ELAISR20J002830-424456B 5.400 
ELAISR20J002830-424456T 

ELAISR20J002830-425102 7.071 

ELAISR20J002830-425349 0.500 

ELAISR20J002835-424822 0.416 

ELAISR20J002837-424528A 1.179 

ELAISR20J002837-424528B 2.510 
ELAISR20J002837-424528T 

ELAISR20J002840-425746 0.455 

ELAISR20J002846-425132 6.822 

ELAISR20J002850-430155 1.025 

ELAISR20J002854-431218 3.567 

ELAISR20J002856-430556 0.383 

ELAISR20J002909-432922 1.111 

ELAISR20J002912-431434 0.841 

ELAISR20J002916-432638 4.683 

ELAISR20J002916-430335 0.516 

ELAISR20J002918-432422 1.081 

ELAISR20J002922-425544 27.700 

ELAISR20J002928-431614 2.893 

ELAISR20J002931-431141 0.611 

ELAISR20J002932-434227 1.619 

ELAISR20J002933-431516 1.043 

ELAISR20J002933-434314 0.403 

ELAISR20-J002934-430811 2.491 

ELAISR20J002938-434447 1.032 

ELAISR20J002938-432946 0.644 

ELAISR20J002945-424142 0.429 

ELAISR20J002946-432149 7.324 

ELAISR20J002947-431554 24.110 

ELAISR20J002951-440551A 2.563 

ELAISR20J002951-440551B 4.419 
ELAISR20J002951-440551T 
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Figure 7. Contour images of the 31 radio sources classified as double or multiple in SI. The contour levels are at 3, 4.5, 6, 7.5, 10, 
20, 40, 60, 80, 160, 320, 640 times the local rms value for all the sources. 



as an example. The full c atalogue will be available from 
http:/ /athena.ph. ic.ac.uk/. 



4.3 Errors in the source parameters 

The formal relative errors determined by a Gaussian fit are 
generally smaller than the true uncertainties of the source 
parameters. Gaussian random noise often dominates the er- 
rors in the data (Condon 1997). Thus, we used the Condon 
(1997) error propagation equations to estimate the true er- 
rors on fluxes, axes and position angle: 




where S p and St are the peak and the total fluxes, 8m and 
9m the fitted FWHMs of the major and minor axes, PA is 
the position angle of the major axis (the as are the relative 
errors) and p is the signal-to-noise ratio, given by 













1 + 






1 + 


( 0jL Y 










\6 M ) _ 






\e m ) 



where a mav is the noise variance of the image and On is 
the FWHM of the Gaussian correlation length of the im- 
age noise (~FWHM of the synthesized beam). The expo- 
nents are om = 5/2 and a m — 1/2 for calculating cm, 
olm = 1/2 an d a m = 5/2 for calculating a m and op a and 
om = a m = 3/2 for the flux densities. These two equa- 
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Figure 7 - continued 



tions are the master equations for estimating the variance 
in the parameters derived from a two-dimensional Gaussian 
fit. The projection of the major and minor axis errors onto 
the right ascension and declination axes produces the total 
rms position errors given by Condon et al. (1998) 

o~l= e 2 a + o~l a sin 2 (PA) + o- 2 yQ cos 2 (PA) (4) 

°l = el + o~l a cos 2 (PA) + a 2 yo sin 2 (PA) (5) 

where (e a ,£s) are the "calibration" errors, while a XQ and 
a yo are 8 2 M /{Aln2)p 2 and 9^/ (Aln2)p 2 respectively. 

Calibration errors cannot be determined from the sur- 
vey image alone. They can be determined from comparison 
with accurate positions of sources strong enough that the 
noise terms in equation W and H are much smaller than the 



calibration terms. Because SI area is covered by no other 1.4- 
GHz radio catalogue suitable for estimating the calibration 
errors, we have used our data to derive the mosaic errors. 
This internal check gives only a lower limit on the source 
position errors. However, Hopkins et al. (1998) have used 
Monte Carlo methods to explore source-fitting uncertain- 
ties in ATCA mosaiced images with statistical properties 
similar to those of the images discussed here. Their study 
implies that the positional errors are rarely as large as 2 
arcsec, and more typically lie below 1 arcsec. 

The mosaic map contains several overlapping regions 
where sources are detected in two different pointings. We 
have run IMSAD on each single pointing map (after hav- 
ing corrected for the primary beam shape), considering the 
sources detected as an independent data set. Then we used 
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Figure 8. Position errors for 36 bright single sources (S > 10 
mjy) common to two individual pointing images composing the 
mosaic image. In the top panel the semi-axes of the 90% confi- 
dence ellipse shown are (2 In 10) 1 / 2 times the rms errors (e a , Eg) 
= (0".138, 0".323). 



Figure 9. Position errors for 236 weak single sources (fainter 
than 10 mjy) in the overlapping regions, in units of the com- 
bined position uncertainties (cr a , ag). In the top panel the 90% 
confidence error circle is plotted. The smooth curves in the central 
and lower panels represent the expected Gaussian of zero mean 
and unit variance. 
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Figure 10. Rms position uncertainties a a and as for all the 
single sources of peak flux S p in our catalogue. 



the sources in common between two different pointing maps 
to estimate the positional errors. First we used 36 sources 
brighter than 10 mjy to calculate the mean image offsets 
and the mosaic uncertainties. Their offsets Aa and AS are 
shown in figure ^. The mean offsets of our mosaic map are 
< Aa >= -0."027 ± 0."029 and < AS >= 0."031 ± 0."059. 
Due to the negligible values of these image offsets, the source 
positions in the catalogue have not been corrected for them. 
The offset distributions (reported in figure ^) indicate the 
mosaic position errors of our data, which are e a = 0."138 in 
right ascension and es = 0."323 in declination. 

Positional uncertainties for all sources have been calcu- 
lated according to equations ^ and ^, adopting the mosaic 
errors derived using the 36 sources brighter than 10 mjy. 
To verify that these uncertainties are realistic, especially at 
low flux densities where the noise and confusion components 
dominate the errors, we used the 236 sources fainter than 10 
mjy in the overlapping regions and found their position off- 
sets. Figure [)] shows these offsets (Aa, AS) normalized by 
the combined uncertainties (a a , as), where a 2 , — a\ x + cr„ 2 
and a s = a si + a S2 with a a i, a a 2, ersi and as2 being the 
errors on the single measurements from equations @and§). 
As shown in figure ^| where the expected Gaussians are plot- 
ted over our normalized offset distributions, the normalized 
offset have nearly zero mean and unit rms scatter, verifying 
that our catalogue uncertainties are accurate also for weak 
sources. Many of these sources are slightly extended, so our 
catalogue positional errors for weak sources include possible 
offsets between the source centroids and cores. In extended 
sources, there might be an intrinsic offset between the fitted 
centroid position and the real radio core position. 

The rms position uncertainties (a a , as) of all the 
sources in our catalogue are plotted as functions of peak 
flux in figure hoi 



Figure 11. The measured ratio of integrated to peak flux as a 
function of signal-to-noise ratio for the SI ATCA survey. The 
upper line defines the upper envelope of the St/Sp distribution 
containing the sources that we have considered unresolved. 

5 SOURCE COUNTS 

The sample of 581 sources with S p > 5a has been used to 
construct the source counts distribution. Complex sources 
or sources with multiple components have been treated as 
a single radio source. Every source was weighted for the 
reciprocal of its detectability area (figure ^|) , defined as the 
area over which the source could have been seen above the 
adopted threshold of 5a (Katgert et al. 1973). 

An estimate of the extension of a source can be obtained 
from the ratio of the integrated flux density to the peak flux 
density St/S p . Figure llll shows that the ratio has a skew 
distribution, with a tail towards high flux density ratios, es- 
pecially when the signal-to-noise ratio is low. To establish 
a criterion for extendedness, we have determined the upper 
envelope of the distribution of St/S p containing the unre- 
solved sources. The solid lines drawn in figure represent 
the best-fit curves to the upper and lower envelopes of the 
band containing all the sources considered to be unresolved. 
For these sources we have adopted the peak flux in comput- 
ing the source counts, while for all the others, lying above 
the upper envelope, we have adopted the total flux. 

In Table |^the 1.4 GHz source counts are presented. The 
columns give the adopted flux density intervals, the average 
flux density in each interval computed as the geometric mean 
of the two flux limits, the observed number of sources in each 
flux interval, the differential source density (in sr _1 Jy 1 ), 
the normalized differential counts nS 25 (in sr" 1 Jy 1,5 ) with 
estimated errors (as n 1 ^ 2 S 2 ' 5 ) and the integral counts (in 
sr" 1 ). 

The 1.4 GHz differential source counts of the ATCA 
data, normalized to those expected in a Euclidean geometry 
by dividing by S~ ' J , are shown in figure |l| (filled stars). 
For comparison, source counts from other surveys are also 
plotted. The solid line represents the model of Windhorst, 
Mathis & Neuschaefer (1990) obtained by fitting counts from 
24 different 1.4 GHz surveys, while the open circles are the 
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Table 2. The 1.4 GHz ATCA Radio Source Counts 
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Figure 12. The 1.4 GHz normalized differential source counts 
for the ELAIS ATCA data. The abscissa gives the flux density 
(mjy) and the ordinate gives the differential number of sources 
normalized by 5 2 s (sr - 1 Jy 1 ' 5 ). The solid line represents the 
model of Windhorst, Mathis & Neuschaefer (1990) obtained by 
fitting counts from 24 different 1.4 GHz surveys. The filled stars 
are the counts obtained from our ATCA data in the ELAIS south- 
ern region, while the open circles are the total radio counts in the 
ELAIS regions obtained by combining the SI data with the VLA 
data in the northern ELAIS regions (Ciliegi et al. 



total counts obtained by combining the SI data with the 
VLA data in the northern ELAIS regions (Ciliegi et al. 
1998). Our counts agree with those obtained by previous 
surveys and confirm the upturn observed by several authors 
below about lmjy, considered as the characteristic feature 
of the 'sub-mjy population'. 

A maximum likelihood fit to both our ATCA and 
ATCA+VLA 1.4 GHz counts with two power laws: 



dN 
~dS 



if s > s b 
if s < s b 



gives the following parameters: Qi = 
0.27, S b =(0.72±0.23) mjy and oti 



1.73±0.11, Q 2 
= 1.74±0.06, a 2 



(6) 

3.04± 
: 2.48± 



0.11, Sb =(0.53±0.13) mjy, respectively. Although the errors 
for the ATCA data only are relatively large, in both cases 
our best fit parameters suggest that the re-steepening of the 
integral counts toward an Euclidean slope starts just below 
~1 mjy, in agreement with the recent results of Gruppioni et 
al. (1997) and Ciliegi et al. (1998). By contrast, in previous 
work (Windhorst et al. 1985, 1990), a best-fit value of 5 mjy 
had been derived for the flux density at which the change in 
slope of the source counts occurs. 



6 CONCLUSIONS AND FUTURE PLANS 

This paper presents a 1.4 GHz survey obtained with the 
ATCA of the European Large Area ISO Survey (ELAIS) re- 
gion SI located in the southern celestial hemisphere. A mo- 
saic of forty-nine separate observations with different point- 
ing positions has provided an image with low and uniform 
rms noise over the whole SI area (~4 deg 2 ). The lowest 5a 
flux density reached by our observation is 0.2 mjy, while the 
bulk of the imaged area has a 5a flux density limit of 0.4 
mjy. 

The observations provide a large, complete sample of 
581 sources (with peak flux greater than 5a), the majority 
of which have flux densities in the sub-mjy range. We have 
constructed the differential source counts over the flux range 
0.2-200 mjy. Our counts agree with the results from other 
deep 1.4 GHz surveys (Windhorst et al. 1990; Gruppioni et 
al. 1997; Ciliegi et al. 1998) and confirm the change in slope 
observed below ~1 mjy. 

This survey of the southern ELAIS field, SI, will be 
combined with the VLA survey of the northern fields, Nl, 
N2 and N3 (Ciliegi et al. 1998), and deep multi- waveband 
data (optical and near-, mid- and far-infrared) in the same 
area, to investigate the nature of the sub-mjy population. 
Moreover, due to the accurate radio positions, the radio sam- 
ple will play a crucial role in the optical identification phase 
of the ISO program. 

Spectroscopic observations for about 350 radio/ISO 
sources in SI have recently been obtained (end of Septem- 
ber 1998) with the Anglo Australian Telescope (AAT) Two 
Degrees Field Spectrograph (2dF). These observations will 
provide crucial information about the spectroscopic nature 
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and the redshift distribution of our objects, and allow esti- 
mation of the obscured star formation rate with two highly 
complementary samples, selected in the infrared and radio 
within the same volume of Universe. In fact, since the ra- 
dio luminosity traces the supernovae associated with star 
formation regions and is not affected by dust obscuration, 
independent estimates of the star formation rate for sources 
with reliable distance indications can also be calculated on 
the basis of radio flux density measurements (see Oliver et 
al. 1998). 
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